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<^|- ' Abstract 

on : 

^ ! The nuclear wobbling motion is studied from a microscopic viewpoint. It is shown 

that the expressions not only of the excitation energy but also of the electromagnetic 
transition rate in the microscopic RPA framework can be cast into the corresponding 
forms of the macroscopic rotor model. Criteria to identify the rotational band associated 

^ ■ with the wobbling motion are given, based on which examples of realistic calculations are 
investigated and some theoretical predictions are presented. 

> : 
^ : 

t— i . SI Introduction 

O; 

The recent advent of new generation crystal ball detectors has been opening a great 
ON ' possibility to explore a new area of the high-spin physics. There are many interesting 
subjects which are waiting to be studied. Among them, we would like to concentrate, in 
this paper, upon the nuclear wobbling motion,^ which is one of the "exotic" rotational 

2 . 

motions in the sense that the axis of rotation does not coincide with any of the inertia 

> : 

r^i ■ axes of deformation. 

The nuclear wobbling motion has been considered by analogy with the spinning motions 
of asymmetric top (classical rigid-body), where perturbations are superimposed on the 
main rotation around one of the principal axes with the largest moment of inertia. When 
quantized, the energy spectra in the energy versus angular momentum plane, which are 
nothing but those of the well-known macroscopic triaxial rotor model (Davydov model), are 
classified into two groups of rotational bands, i.e., the "horizontal" and "vertical" sequences 
corresponding asymptotically to the Regge trajectories associated with the largest and the 
smallest moment of inertia. In the high-spin limit, 1 the physical meaning of those two 
sequences appears to be more transparent by introducing an elementary excitation of the 
"wobbling phonon" mode. The horizontal sequences parallel to the yrast line are rotational 
bands in which zero, one, and two, etc, wobbling phonons are excited in each intrinsic state, 
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while the vertical sequences starting from each yrast state consist of the multiple wobbling 
phonon bands. In fact, the E2 transition rates in the former (|A/| = 2) are larger than 

1 2) 

that in the latter (|A/| = 1) by an order of 1/7 in the high-spin limit. ' 

It is very interesting to ask whether such an "exotic" rotational motion is realized 
as a collective motion in atomic nuclei, because it directly reflects the three-dimensional 
nature of rotational motions. It should, however, be noticed that the concept of the three- 
dimensional rotation is meaningful only when the motion of angular momentum vector is 
traced in the "intrinsic" or the "body-fixed" frame, where the macroscopic rotor model is 
formulated. The very definition, however, of such a coordinate frame is highly nontrivial 
from the microscopic point of view in the general framework of the nuclear many-body 
theory. 3 6) 

Detailed experimental investigation is only possible in the discrete line spectroscopy so 
that the one or two phonons excited horizontal sequences are the most promising targets 
for the study. For such rotational bands near the yrast line, the small amplitude approxi- 
mation to the wobbling mode may be allowed. Then the fully microscopic formulation is 
possible 2 ' 7 10) in terms of the random phase approximation (RPA) on top of the cranked 
mean- field theory, and the analogy to the macroscopic rotor model becomes transparent. 
Thus we mainly concentrate, in this paper, on the yrast and the first excited wobbling 
bands and investigate the possible consequences predicted from the RPA theory in real- 
istic nuclei. In the course of the investigation it will be shown that not only the energy 
spectra but also the interband (AI = ±1) electromagnetic transition probabilities can also 
be expressed in the same way as in the macroscopic rotor model in terms of the micro- 
scopically defined "effective" moments of inertia, which are introduced for the wobbling 
eigen frequency. 

The paper is organized as follows: the E2 transitions in the RPA theory is reviewed 
and is applied to the well-known 7- vibrational bands in §2. Some basic ingredients in 
the macroscopic rotor model, especially the expressions for the electromagnetic transition 
rates, are reviewed in §3 for the completeness. The microscopic formalism is presented 
in §4, while some examples of realistic calculations are studied in §5. Main results are 
summarized in §6. Preliminary results of some part of the present work has already been 
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reported in Refs.. 11 13) 

§2 El Transitions for the 7-vibrational band 

The RPA theory on top of the cranked mean-field approach, which is suitable for the 
high-spin states, was first developed in Ref. 14 and applied to the high-spin (3- and 7- 
vibrational bands. lo l6) The extension to the odd nuclei, with special attention to how the 

17) 

electromagnetic transition rates should be calculated, has been done in Ref.. 

The RPA treatment of the vibrational excitations 14 ' is based on the boson expansion 
theory; the lowest order vacuum states, on which the RPA modes are excited, are described 
by the static mean-field theory uniformly rotating around one of the inertia axes (one 
dimensional cranking states). Combined with the 1/1 expansion technique, the matrix 
element of the electromagnetic transition with multipolarity A is expressed in the simple 
form; 

M(i -> /; AJ) « (/|<2am=a/K>> (2-1) 

where AI = If — Ii and the superscript (E) means the electric part of the transition 
operator Q\^. It should be stressed that the components of Q^J are defined with respect 
to the rotation (cranking) axis. 

In the formula above only the lowest order in 1// is retained. Although we are mainly 
concerned with the high-spin limit and consider only in the order of eq.(2.1) in the fol- 
lowing sections, it is worth while mentioning that the formula is applicable at low-spins 
if the the geometry of the angular momentum vector is properly taken into account. A 

18) 

good example is the Ml transition in odd nuclei at relatively low angular momenta. It 
is, however, noticed that the idea is more general and is based on the observation that 
the most important part of the spin-dependence, which comes from the dynamics not 
from the kinematics of the angular momentum algebra (the Clebsch-Gordan coefficients), 
is contained in the right hand side through the spin-dependent change of cranked wave 
functions. In order to see this is really the case and to show the reliability of the formula 
(2.1), we compare, in the remaining part of this section, the results of the RPA calculations 
for the low-spin 7-vibrational band 16) in even-even nuclei with the experimental data. 
The quantities to be investigated are the E2 transitions between the ground state band 
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and the 7-vibrational band. Therefore, the initial and the final state in eq.(2.1) are 

\i)=Jt^\u vot ), |/) = k ot ), (2.2) 

where |u; ro t) is cranked mean-field approximation of the yrast states with rotational fre- 
quency a; ro t, and Xy ^ is the creation operator of the RPA eigen mode at the corresponding 
frequency, which smoothly continues to the 7-vibrational solution at zero-frequency. The 
superscript (±) 16) denotes the signature quantum number carried by the RPA mode, so 
that, for example, the (+)-band represents the even-spin member and the (— )-band the 
odd-spin if the vacuum |u> ro t) has zero signature. There are five kinds of transitions asso- 
ciated with AI = 0, ±1, ±2, and the S(-E'2)-values of these transitions are calculated by 
means of the electric (i.e. proton) part of the operators, 17 ' 

Q20 = lQ { o +) + ^-Qi + \ (2.3a) 

Q 2 ±i = ^(Qi" ) ±Q<- ) ), (2.36) 

Q2±2 = -^(^Q { +) + IQ ( 2 +) ±Q[ +) ), (2.3C) 

where Q { ^ ] (K = 0,1,2) are the signature coupled quadrupole operators with the z-axis 
as the quantization axis, while the rotation axis is chosen to be the usual. 
Among the five transition amplitudes between the 7 and the ground state band, 

tWP)} = U\Q { k ] \i) = MIQP^^M = ([Q ( k\x^]), (2.4) 

only the K = 2 components are non-zero at u rot = from the selection rule because the 
ground state under consideration is axially symmetric. Moreover, the electric part of them 
are of the form 

t[(Q { K } ) {E) ] = Ql ($K2 + <W 7 c^ot) + OK 2 ot ), (2.5) 

in the low frequency limit. As it is mentioned we have to include the effect of Clebsch- 
Gordan coefficients at the low-spin. In this case (ij22 — 2|i/0) is necessary, which has the 
asymptotic values in the large / limit; 

' y/6/4 AI = 0, 
(/ i 22-2|/ / 0>| w <J 1/2 A/ = ±l, (2.6) 

1/4 AI = ±2. 
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Combining with eqs.(2.3) and (2.6) we have the expression of the B(E2) = \M(i — > /)| 2 , 
which is valid in the low-spin limit, 

B{E2)? 



"y—>g 



|v / 2Q 7 | 2 (l + a 7 A/u; rot + 0(a; r 2 ot )) , (2.7) 



(/<22 - 2\i f oy 

where eq.(2.5) is inserted. Using again the asymptotic relation and the fact that the 
rotational frequency can be related to the angular momentum through the moment of 
inertia J ', 

If (I f + 1) - h(Ii + 1) « 21 AI, hu rot « 7/ J", (2.8) 

where 7 = (7/ + Ii)/2, we finally find 

i 

U K a ( : + 9 lW/ + 1} " m + (2 ' 9a) 

in the first order, which is nothing but the formula discussed as the generalized intensity 
relation in Ref., 19> where 

Q t = |v / 2Q 7 |, q = a J /2hJ. (2.9b) 

It should be emphasized that the parameters Q t and q are now calculated microscopically 
by means of the RPA theory at the finite rotational frequency. 

An example of the RPA amplitudes of the electric quadrupole transition operators, 
eq.(2.4), as functions of the rotational frequency is shown in Fig.l for a typical well- 
deformed nuclei, 164 Er. The procedure of calculation is the same as Refs., 1617> except that 
the difference of the oscillator frequency between neutrons and protons in the Nilsson po- 

20 21) 

tential and of the oscillator length in the quadrupole residual interaction ' are properly 

*) 

treated. Namely, the residual interactions of the monopole pairing and the quadrupole 
interactions are used, the strengths of which are determined so as to reproduce the even- 
odd mass difference, and the excitation energies of the /?, 7 -vibrations (K = 0, 2) and the 



*) In Ref. 16> the RPA amplitudes of the 7-vibration for the mass quadrupole oper- 
ators were shown. For the collective solutions the proper treatment of neutron and 
proton oscillator lengths makes the transition amplitudes of the electric operators 
approximately about Z/A of those of the mass operators, just as in the case of the 
static quadrupole moments, i.e. the mean values of the mass and electric quadrupole 
operators, see eq.(4.32) and discussions at the end of §4. 



zero-frequency (Nambu-Goldstone) mode (K = 1), at u TOt = 0, respectively. After fixing 
the force strengths at uj rot = 0, the pairing selfconsistent calculations are performed as func- 
tions of the rotational frequency, but the deformation parameters are fixed at the values 

22) 

deduced from the experiments, for simplicity. Increasing the rotational frequency, the 
quasiparticle alignments generally occur. In order to identify individual rotational bands 
with the same internal structure, are used the diabatic quasiparticle orbits, 16) specifically 

17) 

the diabatic basis constructed by the method of the w rot -expansion up to the third order. 
The model space of the Nilsson orbits are chosen as A^sc = 4 — 6 for neutrons and = 3 — 5 
for protons, which reproduces the transition amplitude very well without using any kind 
of the effective charge (see below). 

From the figure it is seen that the amplitudes follow the low frequency behaviours of 
eq.(2.5) and the parameters, <2 7 and a 7 , are easily extracted. The values thus obtained for 
some arbitrarily chosen rare earth nuclei are summarized in Table 1, where the parameters 
in eq.(2.9), Qt and q, are also included and compared with the available experimental data. 
Note that the moments of inertia for the ground state band and the 7-vibrational band are 
the same in the low-spin limit within the RPA theory. It is not the case in experimental 
data, however. This is because the difference of the moments of inertia results from the 
higher order couplings between the vibrational and rotational motions, which are not taken 
into account in the RPA. Therefore two values of q which are obtained by using two choices 
of experimentally determined J = J g and J = J 1 are included in the table. 

The absolute values of the amplitudes are well reproduced in our calculations. How- 
ever, the agreement might not be taken so seriously because the resultant RPA amplitude 
depends on the size of the adopted model space. It is known that the calculation with full 
model space usually overestimates the experimental values in the simple monopole pairing 
plus the quadrupole interaction model. We will not discuss this point further in this paper. 

More important is that not only the sign but also absolute values of the q parameter 
are well accounted for in the present RPA theory, which is calculated from the ratio of the 
transition amplitudes and therefore more reliable than the amplitudes themselves. This 
clearly shows the rotation-induced change of the microscopic structure of the vibrational 
motion, i.e. the K-mixing of the transition amplitudes, see eq.(2.5), are correctly described 
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in the RPA theory. In Fig.l the appreciable deviations from the lowest order relations, 
the order of O{u% ot ), are predicted in the region of fuv rot > 0.15 MeV, for example, the 
amount of the reduction of the I 1 — > (I — l) g transition rates are larger than that of 
the 7 7 — > (I — 2) g transitions (see eq.(2.3)). We should, however, be careful to draw a 
definite conclusion: at these moderate spins other higher order effects in 1/1 neglected in, 
e.g., eqs.(2.6), (2.8), are of the same order, and then all the higher order terms should be 
consistently calculated, which is out of the scope of the present calculations. 

§3 Macroscopic Rotor Model 

The nuclear wobbling motion was originally predicted as a collective motion in the 
macroscopic rotor model. 1 ' In order to see the characteristic features and to compare with 
the microscopic model, the energy spectra and the E2 and Ml transition rates in the 
high-spin limit are summarized in this section. 

The hamiltonian of the rotor model is written in terms of the three moments of inertia 
around the principal axes, J x , J y and J z : 

/2 /2 j2 

Hr = — y + Yf + ipfi 

where the angular momentum operators here are components with respect to the body-fixed 
frame. The energy spectrum of the above hamiltonian is well known: it is specified by two 
quantum numbers and they are given in the large / limit explicitly by 

Er(I, n w ) = ^± ^ + hu w (I) {n w + \), (3.2) 
where J x is assumed to be the largest and the main rotation occurs around the x-axis. 

1 23) 

The wobbling frequency is determined by the well-known formula, ' 

hu w = I jWyW z = hu TOt l ^-Wx- M, (3.3a) 

V UyJz 

where 

fko Iot = I/J x , (3.36) 
corresponds to the rotational frequency around the main rotation axis and the quantities, 

W y EE 1/J Z - 1/J X , W z EE 1/J y - 1/J X , (3.3C) 



are introduced. The integer n w = 0,1,2,... in eq.(3.2) is the wobbling phonon number 
excited on the yrast states. Then the horizontal and vertical sequences mentioned in §1 
are precisely 

with n w specifying the yrast, yrare,... bands, and 

E \Z\=iW = Er{ -^ n ™ = I - K )> K = K U K 2 , K 3 , (3.46) 

with K specifying the band head spin. Note that the number of phonons are changed in 
the vertical sequence, while it is unchanged in the horizontal one. These two classifications 
correspond to the band structure connected by the E2 transitions with AI = ±2 and 
AI = ±1, respectively, with the transition energies, in the lowest order in 1/1, 

^(hor) = T 2hu rot , (AI = ±2). (3.5a) 

4 ver) = hu w T huj rot , (AI = ±1). (3.56) 

The E2 transition operator in the rotor model is derived by the basic assumption that 
the quadrupole tensor of the rotor is diagonal in the body-fixed frame; thus 

<> = e(f y{ ±±^D$ + + <-,)}, (3-6) 

and then the transition rates is given, again within the large I approximation, by 

B(E2)lT2 ±2 « (ej) 21 - K " «,) 2 , (3.7a) 

Here the deformation parameters (a y ,a z ) are introduced through the static quadrupole 
moments for later convenience; 

z 2 )p(r)dr = ( 1 -Q^-^Q^) : 

(3.8) 

y 2 )p(r)dr=(Q ( + ) ), 
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in obvious notations. These parameters, (a y , a z ), can be related to the usual " (e 2 , 7)" -like 

*) 

parametrization (€3,7*) which are defined by 

J x 2 k p(r)dr = ^R 2 (l + ^e* cos (7* + y fc)) , (k = 1, 2, 3 = x, y, z) (3.9) 

through 

/ 5 47r / 5 

a v = Y ^ e 2 sin (7* + y), ^ = _ Y^ e 2 sm 7*, (3-10) 

so that they represent the static triaxiality around the corresponding y, z-axes. In the 
same way (a y — a z ) is proportional to the static moments around the main rotation (x) 
axis and then eq.(3.7a) simply denotes the usual stretched E2 transitions with respect 
to the x-axis. If we consider the lowest excited band, the 1-phonon wobbling band, then 
n w = 1 in eq.(3.7b). 

In the same framework Ml transitions can be considered if an appropriate magnetic 

24) 

property of the rotor is assumed. According to Ref., we assume that each component of 
the magnetic moment vector in the body-fixed frame is proportional to that of the angular 
momentum, 

mi = gili (i = x,y,z), (3.11) 

where gi (i = x, y, z) denotes the g-factor with respect to the i-th axis in the body-fixed 
frame. Then the Ml transition operator is derived as^ 

-liay-g^i-iD^ + i^D^)} 

where hn is the nuclear magneton and J^ ab ^ means the angular momentum operator in 
the laboratory frame so that the first term in the curly bracket does not contribute to the 

transitions. The B(M1) is evaluated as in the same way as B(E2); 

1 1 

5(M1)E ±1 - {^ 2 N )nJ- A (( 9y - 9.) (^) \{ 9z - g x ) (^) 'J. (3.13) 



*) These parameters coincide with the original ones only within the first order in 
£2 because the original ones are defined with respect to the anisotropic harmonic 
oscillator frequencies. The Lund convention for the sign of triaxiality parameter 7 is 
used throughout in this paper. 

f) The Ml operator used in Ref. 24 is different from eq.(3.12), which, we think, is 
incorrect. The resultant B(M1) is also different. 

9 



Note that B{E2y^l ±1 oc 1/1 and S(M1)X2 ±1 oc / while B(£2)™ ±2 does not 
depend on / (up to the leading order in 1//), if all the parameters, the quadrupole defor- 
mations, a y and a z , the moments of inertia, J x , J y and J z , and the g- factors, g x , g y and 

are constants against I. Actually, all these parameters are not independent and, for 
example, the quadrupole deformations completely determine the moments of inertia in the 
original rigid-body model. However, if one wants to apply the model to realistic nuclei, one 
should consider that all the parameters depends smoothly on the angular momentum, and 
the moments of inertia do not necessarily take the rigid-body values. As will be discussed 
in the following sections, the microscopically calculated J x , J y and J z change as functions 
of the rotational frequency and the relationship between them are far from that of the 
rigid-body f j) 

It is worth mentioning that the axially symmetric limit should be taken with great 
care. There are two kind of limits: the "collective rotation" limit, where 7 = 0°,— 60°, 
or equivalently either a z = or a y = 0, and the "non-collective rotation" limit where 
7 = 60°, —120°, or a y = a z . In the former limit, it is easy to see that there are no definite 
limiting expressions for the wobbling energy and the AI = ±1 B(E2) and B(M1), as 
long as the detailed limiting behaviours of three moment of inertia are not specified. In 
contrast, one can always argue the latter limit, where, of course, J y = J z so that W y = W z . 
Then, from eq.(3.7), the AI = +1 vertical transition as well as the AI = ±2 horizontal 
transitions are prohibited. In fact, in this axially symmetric limit of the "non-collective 
rotation", the horizontal rotational sequences disappear and only the band heads with 
(I x ) = K are physically meaningful as the vacuum states (not the vacuum band). The 
collective rotation in this case, therefore, is the rotation perpendicular to the symmetry 
axis which is at the same time the non-collective rotation axis (x-axis), and the band 
consists of multiple wobbling phonons with excitation energy 

E^)(K) = hu w (K) + hoj Iot » K/J ± , (J ± = J y = J z ), (3.14) 

where eqs.(3.5b) and 3.3) are used, and with the transfer angular momentum AI = +1 
excited on the high-K band head (vacuum) states; namely, 

E \Z\=iW » + !) " K 2 )/2J ± + const., (I > K), (3.15) 
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consistently to the large / ~ K approximation. Note that the rotational frequency hu) rot = 
K/J x disappears in the final expressions as it should be because it is not the collective 
rotational frequency in this case. This kind of rotational band is specifically called as a 
"precession-band" and often observed as a band excited on top of a high-K isomer state 
and has been well studied microscopically in both realistic nuclei 26 ' 27) and in a schematic 
model. 28> 

From the microscopic view point, there is a collective motion even in the limit of axial 
symmetry of the "collective" rotation, 7 = 0° and —60°: that is nothing but the 7-vibration 
with signature a = 1 which has been just considered in the previous section. The reason 
why the limit does not well-behave is that the "body- fixed" frame or the "principal-axis" 
(see the next section) frame of the quadrupole tensor is not well defined in this limit. 
The calculation in the previous section corresponds to the one in the "uniformly-rotating" 
frame, which is always possible to define. It will be shown, explicitly in the next section, 
that the transformation to the "principal-axis" frame is impossible in the "collective" 
rotation limit from the microscopic view point. 

§4 Microscopic RPA Treatment 

Although the wobbling motion is a kind of oscillatory motion of the rotation axis, the 
shape degrees of freedoms are necessary to consider. Especially what is important is the 
non-diagonal parts of the quadrupole tensor and thus we introduce, 

A 

2 V I 7T 



Qy — Qi ^ — oVa-tt s ^X xz ) a ^ 



T (4-1) 



2 V 4tt 

a=l 

which just appear as the AI = ±1 E2 transition operators in eq.(2.3b). Note that only 
the modes with signature a = 1 are relevant in the RPA order so that only these modes 
are considered in the following. In higher order, however, one should also consider the 
remaining non-diagonal component, Q x = Q^K ' j) with all the RPA eigen modes. 

In the RPA treatment, the in-band AI = ±2 E2 transition does not change the RPA 
phonon number and then the matrix elements in the vacuum (or yrast) band and the 
one-phonon wobbling band are the same. In contrast the interband transition from the 
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one-phonon wobbling to the vacuum band reflects the nature of the wobbling motion itself. 
From the argument in §2 (see eqs.(2.1) and (2.3c)) we have 

B(E2)%^ * (Q™f = (e^) 2 \R\a y -a z )\ (4.2a) 

where we have used eq.(3.8) with an understanding that the deformation parameters are 
those associated with the vacuum states, |c^ ro t)? and 

B(E2)^;l±i « l((Qv(n) =F Q z (n)) {E) )\ (4.26) 

where n means that the n-th RPA eigen mode is considered and the RPA amplitudes are 
defined as 

Q k (n) = (n|Q fc |0) RPA = ( [X n , Q k ] ) (k = y, z), (4.3) 

with X n being the annihilation operator of the n-th RPA eigen mode. Note that Q z is 
anti-hermite while Q y is hermite and this is the reason why the sign in eq.(4.2b) is changed 
from that in eq.(2.3b). Here and hereafter (O) for any operator O means the expectation 
value with respect to the cranked state \uJ ro t)- R is clear that the expression of the in-band 
E2 transition formally coincides with those of the horizontal transition in the macroscopic 
rotor model in eq.(3.7a). How about the interband transition? R is one of the main 
purpose of this section to clarify this point. 

The part of operators, Q y and Q z , relevant in the RPA order can be expanded in terms 
of the RPA eigen modes, 



n:all 



Q 2 = ]T(g 2 (n)Xt-h.c.), 

n:all 

Here (n : all) means that the contribution of the Nambu-Goldstone (NG) mode, 



(4.4) 



should be included: 



Q y {n = nG)=--±=2R 2 a y , 



Qz(n = NG) = — L= 2R 2 a z . 
v 21 
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(4.6) 



Namely the contribution of the NG mode corresponds to the "static" deformation while 

those of the normal mode to the "dynamic" (or vibrational) deformations. It should be 

*) 

mentioned that there exists a kind of "sum-rule" which relates both contributions: 

E Qy(n)Qz(n) = —2a y a z , (4.7) 

n#NG 

which can be easily verified from the identity, 

[Q y ,Q z ] = 0. (4.8) 

Although it is not stated explicitly we are working in the so-called " uniformly rotating" 
(UR) frame in the sense that all observables are based on the vacuum state \ui ro t) which is 
the cranked state uniformly rotating around the one of the principal axes of the deformation 
of the mean-field. The wobbling motion naturally appears in the body-fixed frame, or the 
"principal axis" (PA) frame. It is not a trivial matter to define the PA frame in the 
general framework of the many-body problem. One must introduce "gauge conditions" 3 ' 6> 

2 4) 

which are common in the quantum theory with constraints. According to Refs., ' we 
impose the non-diagonal part of the quadrupole tensor should vanish: 

(Q fc ) pA = (k = x,y,z). (4.9) 

The meaning of these condition are apparent. The PA and UR frame picture are related 
through the Euler angles which are now the dynamical variables, for example, 

( J UR = E Afc(e)(J*) pA , (4.10a) 

k=x,y,z 

fe)uR= E A*(G)I^(e)(Q fe ,)p A , (4.106) 

k,l=x,y,z 

where Qij is the non spherically-coupled representation of the quadrupole tensor, = Q^j 
(ijk-cyclic) and Dij(Q) is the rank-1 D function with the additionally introduced Euler 
angles 6. The three conditions, eq.(4.9), in principle, define the microscopic structure of 
the three Euler angles in terms of the complete set of the observables in the UR frame, 

*) Strictly speaking, it might not be called as a sum- rule because each term in the 
summation in eq.(4.7) is not positive definite. 
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where every microscopic quantity is well-defined, and remove the redundancy between 
the microscopic variables in the PA frame and the collective variables 6. Of course it is 
not so simple: for instance, what is the ordering between the Euler angle operators and 
the microscopic quadrupole tensor in eq.(4.10b)? One must work out very carefully in a 
consistent framework. 3) 

Fortunately, it has been shown that the small amplitude approximation makes the 
situation quite simple. 2 ' 6) The Euler angles can be written in terms of the microscopic 
variables and then the transformation from the UR to PA frame can be done explicitly by 
using them within the RPA order, for example, 



UR' 



(^ x )pA ~(^ x ) 

( lJ y)pA = ( lJ y ~ 9fi>2^ ( 5 2 )ur' 



2# 2 a 2 ^UR' (4.11) 



(*^)pA ~( Jz + 2R 2 a ( 5y)uR- 



I 

K y 

Using the UR frame relations eq.(4.4), the angular momentum operators in the PA frame 
are written in terms of the microscopic RPA eigen modes, 



K) PA = -'E o (ig^-H< 



n/NG 

(4.12) 



namely, the NG mode contribution disappears. Apparently the transformation from the 
UR to PA frame is possible only in the case when a y ^ and a z ^ (see eq.(4.21) below 
for more strict conditions) . Now it is easy to check, with the help of the sum-rule relation, 
eq.(4.7), that (a part of) the commutation relation in the PA frame, the sign of which is 
opposite to the one in the UR frame, holds again within the RPA order, 

[(^»)pa»Wpa]=+ j =<Wpa>- ( 4 - 13 ) 

It should be noticed that the transformation from the laboratory frame to the UR frame 
is unitary while that from the UR to PA is not as is clear from this commutation relation. 
The physical reason why the transformation form the UR to PA is non-unitary is apparent: 
the Euler angles in eq.(4.10) are not simple parameters but are now considered to be the 
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dynamical variables canonically corresponding to the collective angular momenta, which 
is a common feature to the theory with constraints. 3 ' 6 ' 

Next, let us consider the dynamics in the PA frame. The time-dependence in the UR 
frame is governed by the hamiltonian, 

H VR = H- u TOt J x . (4.14) 

Since the Euler angles are the dynamical variables describing the wobbling of the rotation 
axis, the hamiltonian in the PA frame are modified, 

Hp a = H Vt x J x VtyJy Vt z J z 

(4.15) 

~-^UR _ ^yJy ~ ^zJz 

where the operators in the right hand side are in the UR frame (we omit ()ur hereafter) 
fifc (k = x,y,z) are angular frequency operators conjugate to the Euler angle operators 
and within the RPA order, 

O x ^w rot >O y , Q z , (4.16) 

reflecting the small amplitude RPA ansatz. The microscopic structure of the angular fre- 
quency operators are determined by the consistency conditions of the gauge conditions, 2 ' 6) 

i^-(Q k ) pA = ([Q fc ,if PA ]) pA = (k = x,y,z). (4.17) 

Again within the RPA order, we obtain 

z n/NG V 2 y/ 

(4.18) 

where the following results of the RPA eigen value problem in the UR frame are used in 
the second equality in each equation, 

#ur ~ huJ nXtX n = hu n XlX n + f ^{J 2 z - (iJ y ) 2 }. (4.19) 

n:all n/NG 

Eqs.(4.12) and (4.18) show that both the angular momentum vector and the angular 
frequency vector precess around the main rotation axis (x-axis) in the PA frame with the 
amplitudes, 

(Jfc) pA (n) = (n|(J fc ) pA |0) RPA , fi fc (n) = (n\Q k \0) KPA (k = y,z). (4.20) 

15 



For this picture to be consistent the ratios of the dynamic to the static deformation should 
be small, i.e., 

r» = Q y (n)/2R 2 a y = {J z ) pA (n)/I ~ 0(1/ VI) « 1, 

(4.21) 

r z (n) ee Q g (n)/2R 2 a z = - (U y ) pA (n)/I ~ 0(l/v^) « 1, 
Now the "effective" moments of inertia for each RPA eigen mode, which are naturally 



introduced through 



are thus written as, 



hjf (n) ee (J k ) {n)/Sl k {n) (k = y, z), (4.22) 



y w„r„(n) +w rot r z (n)' 



uJ n r z (n) + u) rot r y (n) ' 

where eqs.(4.12) and (4.18) and the definition (4.21) are used. Introducing the quantities 
(c.f. eq.(3.3c)), 

W y (n)=l/Jf(n)-1/J x , 

(J x ee I/hw vot ) (4.24) 

W^ra) =1/Jf(n)-1/J x , 

we have 



ftu; n [ r y (n)/r z (n) 1 =i"Wz(n), 

) v (4-25) 

ftw n ( r z (n)/r y (n) ) =IW y (n), 



from which the well-known wobbling energy formula is obtained as in the same way as in 
eq.(3.3a), 



hu n = I^W y (n)W z (n) = hu rot J }eff (n) j-eff (n) • ( 426 ) 

Note that the right hand side of this equation depends on the eigen mode itself so that it 
is only a formal solution of the RPA equation, to which eq.(4.25) is equivalent. It may be 
worth noticing that eigen value eq.(4.25) can be obtained by the Euler equation for the PA 

2 4) 

angular momentum vector, which can be derived from the TDHF variational principle. ' 
In fact, the Euler equation with the PA hamiltonian eq.(4.15), 

^( J 0pA = £ e «*( J i)pA n *> ( 4 - 27 ) 

jk 
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leads, for the component of the n-th RPA eigenmode, 

-u n (iJ y ) .(n) = fl z (n)I -uJvot(Jz)pAn), 

(4-28) 

-u n (J z ) pA (n) = iQ, y (n)I - u) rot (iJ y ) pA (n), 

which reduces to eq.(4.25) first by substituting Q yjZ (n) = (J y , z ) PA (n)/JiJ'^ E (n) (eq.(4.22)) 
and next by {iJ y ) PA (n) = -Ir z (n) and (j 2 ) pA (n) = I r y (n) (eq.(4.21)). 

Since we are considering the case where eigen-energies of all RPA solutions are real 
and positive, the signs of W y (n) and W z (n) are the same and then the basic eq.(4.25) can 
be formally solved for r y (n), r 2 (n), 

r („)=c J-f^M) 1 ' 
1 (W v {n)\\ 

where a n denotes the sign of W y (n) (= the sign of W z (n)) and the quantity c n is the 
amplitude with which the n-th mode contributes to the sum-rule (4.7), 

r y (n)r z (n) = ^ <=► £ a n c 2 n = 1, (4.30) 

so that C 1 for non-collective solutions and (? n 1 for collective solutions. Using the 
definition of r y {n) and r z (n) (eq.(4.21)), eq.(4.29) shows that the transition amplitudes 
Q y (n) and Q z (n) can be expressed in terms of the static deformation parameters and the 
effective moments of inertia. Inserting them into eq.(4.2b), thus, we finally obtain, 

where we assumed 

Q ( k E \n)~ejQ k (n) (k = y,z), (4.32) 

which are approximately satisfied for the (isoscalar) collective RPA solutions. Now the 
analogy with the eq.(3.7b) (with n w = 1) in the macroscopic rotor model is clear: the 
main difference, except to the assumption (4.32), is the fact that there are many solutions 
in the microscopic RPA treatment and the amplitude 7^ 1 and the sign a n can be 
negative. Note that the negative sign a n = — 1 means J y s z {n) > J x or J y s z (n) < 0, which 
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conflicts our basic assumptions that the x-axis is the main rotation axis and the introduced 
three moments of inertia are physically meaningful, and therefore clearly shows that the 
interpretation of the n-th RPA mode as a wobbling motion is not justified in such a case. 
From this observation we use positiveness of the sign a n , i.e. r y (n)r z (n) > 0, and (? n ~ 1 
as criteria to identify the wobbling solutions among the microscopic RPA eigen modes in 
later applications. 

For completeness, the microscopic expression for the Ml rate corresponding to the 
AI = ±1 E2 transitions, eq.(4.2b), is given: 

£(M1)£2£ ±1 « l(i» y (n) ±^(n)) 2 , (4.33) 

with 

/ 3 

Hk(n) = (n|/^|0) R p A = ( [X n ,J—fJi N (9ih + 9 s s k )) ) (k = y,z), (4.34) 
where g\ and g s are usual orbital and spin g-factors. We here are mainly concerned with the 

17) 

even-even nuclei so that only the vibrational contribution are included. In contrast to the 
case of E2 transitions, where is used the same basic requirement as in the rotor model such 
that the quadrupole tensor is diagonal in the PA frame, there is no simple correspondence in 
Ml transitions between the results of the microscopic RPA framework and the macroscopic 
model. This is because the basic assumption of the Ml operator in the rotor model, 
eq.(3.11), is not always justified from the microscopic viewpoint: for instance the effects 
of quasineutron or quasiproton alignments on the Ml transitions cannot be taken into 

account. If we nevertheless take the same assumption as the rotor model, the operators 

*) 

responsible to the Ml transitions are ifx y oc (g y — gx)(iJy)pA, l^z oc (g z — g x ){J z )vA- 
Then using eqs.(4.11), (4.21) and (4.29) a similar expression to the one in the rotor model, 
eq.(3.13), is easily obtained with the same notice as in the E2 transitions, eq.(4.31). 

In the following let us discuss general features of the high-spin RPA eigen modes. The 
properties of the solutions can be different from case to case depending on the microscopic 
structure of vacuum states. Until up to now, we have used the term "collective" somewhat 
ambiguously. It might be necessary to make it more precise here. The collectivity is 

*) Note that only the component of the magnetic moment vector perpendicular to 
the angular momentum in the PA frame, m — g x I, contributes to the transition. 

18 



usually referred, at least, in two different contexts: (1) the transition matrix elements of 
characteristic observables to the solutions, e.g. the E2 transition for the (3, 7-vibrations 
and the wobbling motion under discussions, are large, and (2) the RPA amplitudes of 
phonon operators spread over many two-quasiparticle states. For example, it is the second 
one that is necessary for eq.(4.32) to be valid. Let us call the former the "large-transition" 
collectivity and the latter the " spread-over" collectivity. The shape vibrational motions in 
the ground state regions almost always satisfy the conditions of both collectivities. From 
our experiences, the spread-over collectivity gets weaker in most cases for the RPA solutions 
excited on the vacuum with rotational aligned quasiparticles, e.g. the s-band, at higher 
spins. One of the reasons for this is that the pairing correlations are reduced by the blocking 

29) 

effects of the aligned particles. In spite of the fact that the phonon amplitudes somewhat 
tend to concentrate on few components, the large-transition collectivity can survive even 
at higher-spins, since the energy lowering of the two-quasiparticle states make the RPA 
eigen energy small so that the backward amplitudes remain to be large. 

The discussions above indicates that the Ml amplitude with using the microscopic 
Ml operator, eq.(4.34), is small for the typical spread-over collective RPA solutions. In 
fact a destructive interference occurs for the Ml operator, because amplitudes of this type 
of solutions, which are of isoscalar quadrupole character, have many two-quasiparticle 
components with phases favourable to the E2 operator. Actually the quadrupole shape 
vibrations in non-rotating nuclei have RPA amplitudes which transfer the K quantum 
number by even unit so that the transition matrix elements for the Ml operator strictly 
vanish. This is well established for the 7-vibrational bands at low spins, where the Ml 
transitions are usually small. Therefore the AI = ±1 transitions are E2 dominant in 
such cases. It should, however, be noticed that the spread-over collectivity of the RPA 
solutions are reduced at higher-spins and then the properties of individual two-quasiparticle 
states, whose components are dominant in the RPA amplitudes, manifest themselves as 
an enhancement of the Ml amplitudes. We will see examples of such cases, where the Ml 
transitions are no n- negligible or even dominant, in the next section. 
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§5 Examples of Numerical Calculations 

In the microscopic RPA formalism we have solutions of vibrational character in many 
cases. As it is clarified in the previous sections, however, all of the solutions cannot 
necessarily be interpreted as the wobbling-like motions even though they are collective 
enough in the sense that the E2 transition amplitudes are large (see discussions below 
eq.(4.31)). Namely, in addition to the enough collectivity, the relative sign of the transition 
amplitudes Q y (n) and Q z (n) should be the same as that of the static asymmetries around 
y and z-axes, a y and a z , 

sign of (Q y (n)/Q z (n)) = sign of (a y /a z ). (5.1) 

This relation was mentioned in another form in eq.(4.9) of Ref.. 25) Only in such solutions 
the three moments of inertia associated with the RPA solutions, eq.(4.22), are well defined 
and are consistent to the picture of the wobbling motion. Interestingly enough, the lowest 
collective RPA solutions, if they exist, satisfy the condition in most of the cases we found. 
It might be interesting to note that eq.(5.1) holds also for 7-vibrational excitations in 
nuclei with small selfconsistent 7 deformation: It is a good approximation to set 7 = 
in such nuclei and then the same relative sign as the selfconsistent triaxiality (eq.(5.1)) is 
obtained for the transition amplitudes even with such an artificial setting of 7 = 0. 

More precisely speaking, one must also require the conditions eq.(4.21) in order for 
the adopted small amplitude approximation to be valid. This means that both static 
deformation around y and z-axes should be larger than the dynamic ones. Therefore, one 
might expect the ideal wobbling motion only in the region of equilibrium shape with 7 7^ 0° 
and 7 7^ —60° in the (e2,7) plane; i.e. an appreciable amount of triaxiality is generally 
required. Notice that the wobbling-like solutions can exist in the axially symmetric cases 
with 7 = 60° and = —120°, where a y = a z < and > 0, respectively, and correspond to 
the AI — 1 excitations on a high-K isomer states (the precession-bands) as is discussed in 
the end of §3. 

This apparent asymmetry between the four cases of axially symmetric shapes are based 
on the fact that the three rotation axes are not treated equivalently: the x-axis is the main 
rotation axis and the angular momenta around the other two axes are small compared 
to it, which is generally believed to be valid near the yrast states. Such a treatment 
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especially makes sense in the situation where the quasiparticle alignments occur: then an 
appreciable part of the main rotation are carried by single-particle degrees of freedoms, 
while perpendicular components of the rotation are only of collective nature. Actually, we 
have found the wobbling-like solutions only in the cases where aligned quasiparticles exist 
in the vacuum state from which the RPA modes are excited, which makes the vacuum 
time-reversal broken and induces strong rotational K-mixing effects. 

With the general considerations above in mind we will show some examples of the 
results of realistic calculations and discuss the characteristic properties of the obtained 
wobbling motions. 

§5-1 AI = ±1 interband E2 transitions 

Restricting ourselves to even-even nuclei, the yrast band has naturally the quantum 
numbers, (a, tt) = (0,+) with even-spin values. Then the (one-phonon) wobbling band 
has (a,7r) = (1,+) with odd-spin values. Unfortunately, there are only few cases where 
(a, 7r) = (1,+) bands are observed up to high-spin states. Moreover it is generally difficult 
to distinguish the collective phonon band from non-collective quasiparticle bands, like the 
AC or BD two-neutron aligned bands in the usual nomenclature of quasiparticle orbits, 
only from the energy spectra. Therefore, at least until now, there is no definite evidence 
that the ideal wobbling motions exist in atomic nuclei. 

One of the most important feature of the wobbling motions is that the E2 amplitudes 
Q y (n) and Q z (n) have comparative magnitude, in contrast to the case of the 7- vibrations 
in §2, where \Q y (= t[Q[ ' ) ])| <C \Qz(= — '*])!• This property immediately leads that 
the one of the transitions with either AI = +1 or = —1 is much larger than the other 
depending on the relative phase of Q y (n) and Q z {n), see eq.(4.2b), or the triaxiality, see 
eq.(5.1). This is situation quite analogous to the Ml transitions between the signature 
partner bands in odd nuclei, where the single-particle matrix elements of the operators i\x y 
and \i z between the signature partner bands play a similar role to the two amplitudes Q y 
and Q Z: see Figure 2, although the role played by the triaxiality is not necessarily the same. 
For the ideal case of the wobbling motion where the formula of the macroscopic model are 
valid, the staggering of the E2 transitions between the AI = +1 and —1 combined with 
the excitation energy gives important information on the triaxial deformation and/or the 
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three moments of inertia; for example, 

AB(E2) ]AI \ =1 B(E2) AI= _ 1 - B(E2) AI=+1 _ 8 a y a z 
B(E2) lAI \=2 B(E2) lAIl=2 I(a y -a z ) 2 ' [ ' 1 

just as in a similar kind of analysis for the staggering of the Ml transitions in odd nuclei as 

in Ref.f 0> The relation between the B(E2) A j = ±i and the triaxial equilibrium deformation 

is summarized in Figure 3. 

25) 

Although there is no definite evidence, it is suggested in Ref. that the odd-spin se- 
quence of the so-called "extention of the 7-band" 31) after the g-s band crossing in 182 0s 
might be a candidate of the wobbling motion from the calculations of the same RPA for- 
malism. It is instructive here to show the results for this nucleus because the results of 
the calculations show an ideal feature discussed in the previous sections. This nucleus 
is supposed to be 7-soft and the two quasineutron alignment induce the negative 7 de- 
formation. Although the potential energy surface are rather flat so that it is difficult to 
determine the precise 7 value, the RPA results with 7 = —18° has been shown to have 
desired property as an ideal wobbling motion. 12 ' 2j) It should be stressed that in order to 
obtain the wobbling-like RPA solution the existence of the aligned quasiparticles in the 
vacuum configuration, e.g. the s-band, is indispensable. 

The calculational procedures, the single particle potential, the effective interactions 
and the model space etc, are the same as those in §2. It should, however, be stressed 
that the RPA dispersion equation with the NG modes being explicitly decoupled 16 ' is used 
in this calculations. For the case with 7 7^ 0°, then, the quadrupole force parameters 
disappears for the signature a = 1 sector. Therefore the wobbling-like RPA solution are 
obtained without any ambiguity depending on the choice of the force parameters. Just 
as in Ref., 2j) *' ) all the mean-field parameters are fixed, for simplicity, with the values 

*) In Refs., 12,2o> the transition amplitudes are 10 — 30% larger because a larger model 
space (N osc = 5 — 7 for neutrons and = 4 — 6 for protons) than that adopted in this 
paper has been used. From the results in §2 we hope the calculations in this paper 

25) 

are more realistic. In addition, the results in Ref. were obtained by using the 
RPA equation without the NG-modes explicitly decoupled. The rotational frequency 

dependence of the excitation energy of the wobbling is thus a little bit different form 

12) 

that in Ref., where the same NG-decoupled equation as the present paper was 
used. 
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A„ = 0.74, A w = 0.98 MeV and e 2 = 0.21, 7 = -18°, which are roughly appropriate 
for the s-band of 182 Os. The results of the RPA amplitudes for the lowest RPA solution 
and the corresponding B(E2) values are shown as functions of the rotational frequency in 
Figure 4. Since —60° < 7 < 0°, a y < and a z > 0, and then the relative sign of Q y (n) 
and Q z (n) for the solution is negative (see eq.(5.1)). Thus, the transitions from (J — l) wo b 
to (I) s is much stronger than the one from (/ + l) wo b to (I) s , as is clear in Fig. 4 (see also 
Fig. 3), so that the B(E2) values show remarkable zigzag behaviour. 

Unfortunately these transition rates have not yet been measured. Since the 7-ray 
energy for the (/ — l) wo b — > (I)s transition is much smaller, the 7 decay rate for this 
transition, which has larger B(E2), is hindered. This might be the reason why it is difficult 
to measure the transitions. As clear from the general consideration, if the wobbling type 
RPA solution exists for the positive 7 equilibrium shape, the relative sign of Q y (n) and 
Q z (n) is positive and then both the B(E2) and the 7-ray energy for the (I + l) wo b (-0s 
transition is larger. Therefore, it might be more easy to measure the transitions in this 
case, although we could not find a good example of such calculations. 

The Ml transitions between the wobbling and the s-band in this nucleus is small, 
typically B(M1)(ij, 2 n )/ B(E2)(e 2 b 2 ) pz 0.1, except the highest frequency shown in Fig. 4, 
where Ml amplitudes start to grow rapidly. The smallness of the Ml amplitudes comes 
from the fact that the wobbling solution in this nucleus is an ideal case and keep large 
collectivity up to rather high-spins (see the discussion at the end of §4), although it reduces 
gradually as is seen from the S(i?2)-values in Fig. 4. 

Finally, it should be mentioned that the second RPA solution in this nucleus is also 
considerably collective. 20 ' However, the relative sign of Q y (n) and Q z (n) is positive and 
therefore it is a kind of vibrational mode but is not of wobbling nature. 

$5-2 Effective moments of inertia 

Another important outcome of the microscopic RPA formalism is that the three ef- 
fective moments of inertia can be calculated, which are highly nontrivial from the mi- 
croscopic viewpoint. It should be noticed that the ratios J y s {n)/J x and J z n {n)/J x 
can be extracted from the ratios u n /uj TO t and Q y (n)/Q z (n) through eqs.(4.26) and 
(4.29); the latter ratios are experimentally observable from the energy spectra and the 
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S(£?2)| A j| =1 /S(£?2)| A 7| =2 . It may therefore be interesting to see how the three moments 
of inertia calculated from the wobbling RPA solution behaves as functions of deformation 
parameters, especially the triaxiality 7. Their dependence on the rotational frequency have 
been already studied in Ref. 2 ° for 182 0s, and will be presented for 124 Xe in the next sub- 
section. In the case of the quadrupole residual interaction, a simple formula for Jy S z (n) in 
terms of the 2x2-coupled RPA dispersion equations exists. 16,2j> We have used this formula 
in the following calculations in place of the original definition, eq.(4.22). 

Of course the true equilibrium shape has a definite 7 deformation. Here, however, 
we fixed the other parameters of the mean-field potential and have performed the RPA 
calculations with changing the 7 deformation as a free parameter. Thus, it should be 
considered to be a kind of theoretical simulations how the nucleus behaves if the triaxiality 
is artificially changed. The results for the three sectors in the (€2,7) plane, —120° < 7 < 
—60°, —60° < 7 < 0° and 0° < 7 < 60°, in each of which a representative nucleus, 
176 Hf, 182 0s and 148 Gd, is respectively selected, are shown together in one panel in Figure 
5. Note that the yrast sequences of 176 Hf and 148 Gd are known to be composed of non- 
collective rotations with prolate and oblate shapes, i.e. 7 = —120° and = 60°, respectively. 
Since all the other mean-field parameters are fixed with neglecting the selfconsistency, the 
wobbling-like solutions are not obtained for all values of 7 deformation, and, moreover, the 
solutions can be discontinuous as functions of 7 because of the existence of virtual level 
crossings of the quasiparticle orbits when changing 7. 

In spite of the deficiencies of this relatively simple-minded calculations, it is instructive 
to see the microscopically derived 7 dependence of the three moments of inertia, which is 
neither of irrotational like nor of rigid-body like. In these calculations all the vacuum states, 
on which the RPA mode excited, have aligned quasiparticles, which is essential to obtain 
the wobbling-like solutions as is mentioned above, so that the J x = I /%u rot (kinematical 
moment of inertia) changes very gradually and never vanish. Thus the 7 dependence of 
the three moments cannot be that of irrotaional. On the other hand (k = y, z) are 

dynamical moments of inertia and take as small values as zero. Thus the 7 dependence of 
the three moments cannot be that of rigid-body, either. 

It is very interesting to see how the three moments of inertia, which can be experi- 
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mentally extracted from the observed wobbling motions, if they exist, behave. For this 
purpose, the triaxial deformation should be also determined from the independent exper- 
imental observables, e.g. by the Ml transitions in neighbouring nuclei. 

§5-3 Wobbling motion in 124 Xe 

As pointed out above, one of possible candidates of wobbling motion might be the high- 
spin continuation of the odd-spin members of the so-called 7-band. Actually a scenario of 
character change from the 7-vibrational to the wobbling-like band was first pointed out 
in Ref., and such a trend, though not well developed because of the small triaxiality, 
are suggested for 164 Er in Ref.. 16> The essential point is that a structural change of the 
vacuum caused by the quasiparticle alignment transfers an appreciable amount of K = 1 
quadrupole strength from the Nambu-Goldstone mode to the normal modes. Although it 
is difficult to predict precisely in which nuclei wobbling mode appears because this mode 
is an outcome of subtle interplay between rotation, quadrupole and pairing correlations in 
realistic nuclei, we have at least up to now two typical candidates of wobbling motion in 
observed bands with negative 7 deformation: One is in 182 0s after the (i^iis/2) 2 alignment 

discussed in Ref. 2o> and above, and the other is in 124 Xe 32) after the (z^/iii/2) 2 alignment 

*) 

being discussed in the following. An important characteristic common to both nuclei 
is that the low-spin spectrum is almost rotational but relatively 7-soft while negative 7 
deformation, -30° < 7 < 0° in 182 0s and -60° < 7 < -30° in 124 Xe, is stabilized above 
the first band crossing caused by the high-Q quasiparticles. Here we note that in 182 Os 
the 7-ray which links (the candidate of) the odd-spin members of the 7 band above and 
below the band crossing has not been observed. From this fact, we imagine it might well 
be the case that there might exist wobbling-like bands but have been assigned incorrectly 
as non-collective two-quasiparticle bands in other nuclei because of the lack of the linking 
transitions. A possible mechanism of this lack will be discussed later in §5-4 for the case 
of 126 Ba. 

The low-spin nuclear structure of Xe isotopes has been studied extensively by means 

33) 

of the interacting boson model and known to depend smoothly on the neutron number. 

34) 

A recent study, however, pointed out that above the first band crossing there are some 
*) We thank R. Wyss for informing us of these data. 
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properties which show sudden changes between A = 120 and 122 due to the shape coexis- 
tence; the lighter isotopes remain nearly axially symmetric while the heavier ones become 
triaxial with 7 < 0. In particular, in 120 Xe, a candidate of the even-spin continuation 
of the 7 band above the second band crossing has been found for the first time, whereas 
in 124 Xe, the odd-spin sequence has been known to extend up to higher-spins than the 

32) 

even-spin one and it shows well-developed wobbling character. We mainly concentrate 
on 124 Xe in the following, although similar results can be obtained for other isotopes with 
similar triaxial deformation. 

The result of the RPA calculation for 124 Xe based on the (^/in/2) 2 aligned configuration 
is shown in Figure 6 together with the experimental data seen from a reference which 
makes the routhian of the s-band flat. The method of the calculation are the same as 
in the previous sections but the smaller model space, AT osc = 3 — 5 for both neutrons 
and protons are used in this region of nuclei. This calculation have been done with using 
the deformation parameters £2 = 0.19 and 7 = —45°, which are taken from the Total 

35) 

Routhian Surface (TRS) calculation at huj ro t = 0.3 MeV in Ref.. The hexadecapole 
deformation is neglected for simplicity. Although the meaning of the shape parameters is 

35) 

slightly different because of the different single-particle potential adopted in Ref., we 
confirmed the stability of results of our calculations against possible small change of shape 
parameters. The pairing gaps =1.1 MeV and A„ =0.9 MeV typical for the low spin 
part of the neutron s-bands in this region are used. This RPA calculation reproduces the 
data very well. Note, again, that we do not have any adjustable parameter in the step of 
RPA for the signature a = 1 sector. Calculated effective moments of inertia are shown in 
Figure 7. They vary as functions of the rotational frequency in spite of the fact that the 
nuclear shape is fixed as in the case of 182 0s, 2j) but the dependence is much weaker in 
124 Xe. 

Calculated E2 and Ml transition amplitudes are shown in Figures 8 and 9, respec- 
tively, as functions of the rotational frequency, where the standard values of g\ and gi iree ^ , 36> 
and gi eS) /gi bee) = 0.7 are used in this and the following calculations for the Ml opera- 
tors, eq.(4.34). In the figures, the transition probabilities are also included. Since the 
selfconsistent 7 = —45°, the relative phase between Q y and Q z is negative so that the 
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AI = +1 transitions are stronger as in the same way as in 182 Os. It is, however, noticed 
that the collectivity of the RPA solution, especially the spread-over one, is weaker than 
that in 182 0s, and consequently the Ml transitions is stronger in 124 Xe. Actually, as is 
seen from Figs. 8 and 9, Ml transitions can compete with E2 depending on the transition 
energies. The zigzag behaviour of the E2 and Ml transitions are the same in this case, 
see Figs. 8 and 9. It might be interesting to point out that trends that the E2 transitions 
decrease against the spin while the Ml transitions increase coincide with the prediction of 
the macroscopic rotor model in this case, see eqs.(3.7b) and (3.13). 

We use the quadrupole residual interaction which has two components, K = 1 and 
K = 2, in the signature a = 1 sector. Note that a triaxial deformation mixes K quantum 
number by two unit so that the K = 1 and 2 components of the interaction couple only 
through the rotational motion. It is the interplay of these two components intermediated 
by the rotational coupling that brings about the character change of the lowest collective 
excitation. To look at this mechanism more closely, it is useful to see the second-lowest 
RPA solution together with the lowest one that we have been discussing. The excitation 
energies of the two solutions and the K = 1 and 2 quadrupole transition amplitudes, 
Q y (— t[Qi~^]) and <2 2 (= — t[Q^]), associated with them, calculated at Tiu ro t = 0.3 MeV 
are shown in Figure 10 as functions of 7 from —60° to 0°. The calculation have been 
done using the same mean-field parameters (except 7) as in Fig. 6. The second-lowest 
solution shares an appreciable amount of Q y (n) and Q z (n) strength with the lowest one 
but with the different relative sign between them. The similar situation has already been 
encountered in the case of 182 Os . 2j) In order to interpret these results, here we consider 
the axially symmetric limits, 7 = 0° and —60°, where the symmetry axes are z and y 
axis, respectively. Since the role of y and z axes are interchanged (with the rotation axis 
unchanged) in both cases, the modes with \Q Z \ ^> (<c) \Q y \ in the 7 = 0° on one hand 
corresponds to those with \Q y \ 3> (<C) \ Q Z \ in 7 = —60° on the other hand. We can 
see clearly in Fig. 10 that this correspondence is actually holds and, moreover, that the 
strongest .fT-mixing of the two amplitudes occurs at around 7 = —30°. If we call the mode 
with \ Q Z \ 3> \ Q y \ in the 7 = 0° "7-vibration like", then so is the second- lowest solution in 
124 Xe at the rotational frequency shown. In contrast the lowest solution is "7-vibration 
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like" in 182 0s at hu rot < 0.25 MeV as is shown in Fig. 4. It should, however, be noticed 
that the lowest solution follows the "phase rule", eq.(5.1), which allows us to interpret the 
lowest solution as a wobbling motion in both cases. According to the previous discussions 
the AI = +1 E2 transitions are larger for the lowest solution because (Q y / Q z ) < 0, while 
the AI = —1 E2 transitions are larger for the second-lowest solution. No candidates of 
the second-lowest collective excitation has been observed so far. According to the TRS 
calculation, 35 ' the equilibrium shape of the (i^/i 11( / 2 ) 2 band of 122 Xe is 7 ~ —30°, where 
the excitation energy of the second-lowest solution becomes low and the difference between 
B(E2)ai=±i becomes conspicuous due to the X-mixing. So we think 122 Xe is a more 
promising nucleus for which our theoretical prediction can be tested. 

$5-4 Continuation of 7 -band in 126 Ba 

Next we study the high-spin continuation of the 7 band in 126 Ba. 37) This nucleus is also 
relatively 7-soft at low spins, while the first band crossing is caused by the low-fi (7r/in/ 2 ) 2 
in contrast to the isotone 124 Xe studied above. See Fig. 6 of Ref. 3o> for a summary of the 
data. The alignment of high-j, low-fi quasiparticles drives the nucleus towards 7 > 0. We 

have performed the RPA calculation for excitations on top of the (7r/in/ 2 ) 2 configuration 

37) 

adopting e 2 = 0.24 and 7 = according to the TRS calculation just as in the case of 
124 Xe. On the other hand the pairing strengths are chosen so as to reproduce their band 
crossing frequencies (Fig. 12 of Ref. 37> ) both for proton {Tiuo c ~ 0.35 MeV) and neutron (~ 
0.44 MeV) when the gap selfconsistent calculation has been done. We also confirmed that 
the parameters thus adopted reproduce the observed B(E2 : 2+ — > 0^~) = 1.9 ± 0.2e 2 b 2 at 
the ground state. The calculated pairing gaps for the proton s-band are 1.16 > A„ > 1.03 
MeV and 0.76 > A n > 0.71 MeV for 0.4 < hu rot < 0.5 MeV. Since 7 « 0° in this case we 
cannot use the RPA equation with the NG-mode fully decoupled. We therefore have used 
the equation with the NG-mode partially decoupled 16 ' suitable for 7 = 0°, where the force 
parameter % 2 ^ remains and has been fixed so as to give the correct excitation energy of 
the 7 vibration on the ground state as %uJ 1 (huj rot = 0) = 0.873 MeV. 

The result of the calculated routhian is shown in Figure 11 together with the data 
seen from a reference which makes the routhian of the s-band flat, while the E2 and Ml 
amplitudes are depicted in Figures 12 and 13 with corresponding transition probabilities 
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included. The excitation energy decreases as the rotational frequency increases both in the 
data and calculation; this is contrary to the macroscopic wobbling motion. This result, 
together with the fact that the calculated signature-splitting of energy between a = 
(not shown) and = 1 RPA solutions is very small, favours an interpretation that the band 
under consideration remains to be 7-vibration like rather than evolving to wobbling like 
above the band crossing, as in the case of 164 Er. 16) In fact, the mixing of the K = 1 
component is not so strong, as it is shown in Figure 12, which supports this interpretation. 
The collectivity measured by the sum of the squared backward amplitudes is similar to the 
124 Xe case; ~ 0.6 in both nuclei. Consequently the magnitudes of the interband B(E2) are 
also similar, although the AI = — 1 B(E2) is now larger than the AI = +1 one because 
Qy/Qz > 0, which is characteristic to small positive-7 nuclei, see eq.(5.1). This trend of 
the triaxiality reflects the shape driving effect of the low-fi aligned (7r/iii/2) 2 - The zigzag 
behaviour of the Ml transitions are not pronounced at all in this case. This is because 
the \i z amplitudes are very small. In addition, Ml transitions is decreasing as a function 
of the rotational frequency, which is in contradiction to the prediction of the macroscopic 
formula. 

An interesting feature of the data in this 126 Ba nucleus is that two branching ratios for 
decays from the members of the band under consideration, T(17+, — > 16+)/T(17+, — > 15+ ) 
and T(19+, — > 18+)/T(19^ y — > 17+,), where s and S7 stand for the s-band and the 
a = 1 7-band excited on top of the s-band, respectively, were measured. Using the 
calculated values B(E2 : I S1 -> (I - l) s ) = 0.044 e 2 b 2 = 12 W.u. (c.f. Fig.12) and 
B(E2 : I S1 -> (J - 2) S7 ) ~ B(E2 : (/ - l) s -> (J - 3) s ) = 0.69e 2 b 2 = 184 W.u. at 
huj TOt = 0.4 MeV and measured 7 ray energies, the estimated branching ratio is T(E2 : 
i" S7 — > (I—1) S )/T(E2 : I S1 — > (/— 2) S7 ) ~ 0.05, which is an order of magnitude smaller than 

37) 

the measured one. This apparent contradiction between the data and the calculation 
can be solved by the specific property of the RPA solution in this case such that the Ml 
transition is very strong, as is shown in Fig. 13. Compare the B(M1) values with those of 
124 Xe in Fig.9. 

As is pointed above the collectivity of the calculated solution in this case is more or less 
similar to that in 124 Xe, which is a little bit weaker than that in 182 Os. Actually, looking 
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into the details of the microscopic structure of the RPA phonon, we have found appreciable 
amount of concentrations of the phonon amplitudes to the two-quasiparticle components 
with (7r/i n / 2 ) 2 in 126 Ba and with (vhn/2) 2 in 124 Xe, in contrast to the strong spread-over 
collectivity seen in the case of 182 0s. The concentrations to the different orbits in 126 Ba 
and 124 Xe comes from the fact that the aligned quasiparticles in the vacuum configurations 
are protons in the former and neutrons in the latter, and then clearly explains why the 
Ml transition is an order of magnitude stronger in 126 Ba than in 124 Xe (note also that the 
amount of concentration in the former is a little bit stronger than in the latter). In fact 
the value B(M1 : I S1 -> (I - l) s ) = 0.28 //gr = 0.15 W.u. calculated at %u vot = 0.4 MeV 
(c.f. Fig. 13) explains the observed branching ratio very well. We think, therefore, these 
126 Ba data are the first evidence of the rotationally induced "collective" Ml transition 
between the 7- and yrast bands in the sense that it has much larger Ml transitions than 
E2 transitions though it is an even-even nucleus. 

This Ml matrix element on the other hand gives a possible explanation why the 
linking transition between the S7 and (77 bands in the a = 1 sector was not observed. The 
calculated reduced transition rates give T(M1) = 0.23 x 10 13 s _1 and T(E2) = 0.015 x 
lO 13 *- 1 for 15 S7 -> U a {E 1 = 0.781 MeV) while T(E2) < 0.005 x lO 13 ^ 1 ^ for 15 S7 -> 
13 37 (E 1 = 0.354 MeV). The latter should be multiplied by a hindrance factor due to 
the band crossing in reality. Consequently the intensity of the S7 band flows to the s- 
band rather than the (77 band. The two mechanisms contribute to this results: the small 
transition energy in I SJ — > (J — 2) S7 and the large Ml transitions in I S1 — > (I — l) s . In 
the case of 182 0s with negative 7 deformation, however, this mechanism does not seem to 
apply because i? 7 (9 S7 — > 8 S ) is also small and the Ml transition are hindered. 

§5-5 Effects on properties of odd- A nuclei 

As has been discussed in the previous sections, the character of the lowest lying col- 
lective excitations in the a = 1 sector is 7-vibration like below the first band crossing 
while it sometimes acquires appreciable K = 1 quadrupole strength and therefore becomes 
wobbling like above the band crossing. Rotational bands based on high-j one-quasiparticle 

*) This upper limit is estimated assuming a pure rotational E2 transitions as if no 
band crossing would exist. 
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configurations in odd- A nuclei decouple into a pair of AI = 2 sequences labeled by the dif- 
ferent signature quantum number. The vibrational excitation mode under consideration in 
even-even nuclei affects signature-dependent properties of adjacent odd- A nuclei through 

17) 

the particle-collective coupling. We briefly summarize here the basic consequences of 
the dynamic coupling effects. 

B(E2 : J — > J — 1) connecting these two sequences has been known to show signature 
dependence, i.e., zigzag behavior as a function of spin, which has been observed experimen- 
tally in some one-quasiparticle bands 38 ' 39 ' before the band crossing. The theoretical studies 
has revealed that the signature dependence is brought about not only by the static 40 ' but 
also by the dynamic 41 43> triaxial deformations. These two effects usually contribute in 
opposite sign, and, for example, the effect of the dynamic one is partly canceled by that 
of the static one in 157 Ho, 38> while the latter is dominant in 161 Dy. 39 This dynamic effect 
is nothing but the manifestation of the K = 2 collectivity of the vibrational excitation 
considered in this paper and plays an important role to understand the electromagnetic 
properties of odd- A nuclei. 

Above the band crossing, i.e. in three-quasiparticle bands (e.g. (tt/i 11 / 2 ) 1 (^H3/2) 2 ) m 
odd- A nuclei, the signature splitting of the excitation energy between the two sequences 
often becomes smaller than before the crossing and even becomes inverted in some cases; 
the so-called signature-inversion. This phenomenon has been explained in a similar way to 

44 — 48) 

that for odd-odd nuclei as a result of positive 7 deformation. But this mechanism is 
not appropriate to explain the inversion in the three-quasiparticle bands which has been 
supposed to have negative-7 deformation, for example, in the cases of 165 ' 167 Lu. 9) One 
of the present authors showed that the phase rule, eq.(5.1), results in stronger particle- 
collective coupling in the unfavoured (/ = j + odd) sector than in the favored (I = j + even) 

50) 

sector in negative-7 nuclei and consequently causes signature inversion. In this model 
the effect of the static triaxiality is usually stronger and therefore the explanation in 

44 43 ) 

Refs. consistently survives for positive-7 nuclei. This example clearly shows, again, 
the importance of the collective vibrational modes in the high-spin spectroscopy. 
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§6. Concluding Remarks 

We have studied the nuclear wobbling motion from a microscopic viewpoint with pay- 
ing special attention to the electromagnetic transition properties, especially E2 and Ml, 
and to their implications in relation to the macroscopic Bohr-Mottelson model. Our ba- 
sic stand point is the microscopic RPA theory suitable for high-spin states, which has 
been proposed by Marshalek 214) and developed by the present authors. 1617) In order to 
check the reliability of the microscopic framework, the rotational perturbation to the E2 
transitions between the 7-vibrational and the yrast bands in the low-spin region has been 
examined. The Bohr-Mottelson's generalized intensity relation are naturally comes out 
with the intrinsic transition moments which can now be calculated microscopically. The 
results of the calculations show both the absolute magnitude and the relative phase of 
these moments are well reproduced for some typical rare earth nuclei. 

The strong rotational perturbation at high-spin states causes the structural change of 
the vacuum configuration, i.e. the quasiparticle alignments or the band-crossings. Accord- 
ingly the vibrational excitation modes are naturally expected to change their characters. In 
order to understand the character change, the interpretations of the observable quantities 
are necessary. For this purpose the microscopic wobbling model based on the RPA the- 
ory has been reinvestigated in the light of the macroscopic wobbling model. 1 ' It has been 
shown that the expressions not only of excitation energy but also of the electromagnetic 
transition rates, especially the AI = ±1 E2 transitions, can be cast into the form similar 
to those given in the macroscopic rotor model. Moreover, besides the collectivity of the 
RPA solutions, a criterion to interpret the wobbling motion has been clarified: The relative 
sign of the dynamic quadrupole asymmetries around the axes perpendicular to the main 
rotation axis, i.e. the transition amplitudes Q y (n) and Q z (n), eqs.(4.1) and (4.3), has to 
be the same as that of the static asymmetry around each axis, a y and a z , eq.(3.8). Just as 
in the rotor model, the E2 transitions depend also on the three moments of inertia which 
determine the excitation energy of the wobbling band at the same time. Therefore the 
wobbling motion reflects the three dimensional nature of the nuclear collective rotations. 

The wobbling motion is an unique rotational motion in the sense that the rotation axis 
deviates from the inertia axis of the nuclear body. Recently, similar rotation scheme, called 
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"tilted axis cranking", is proposed and has been applied to the study of the rotational 
bands with no-signature splitting. Although both rotations are of the three dimensional 
nature, they are conceptually different: the wobbling motion is not the stationary motion 
and the angular momentum and the angular frequency vectors draw different trajecto- 
ries in the body-fixed frame, while the tiled cranking is the stationary rotation and then 
the angular momentum and angular frequency vectors are parallel each other with fixed 
direction with respect to the body- fixed frame. In fact, these two vector have a definite re- 
lationship in the wobbling motion, which is nothing but determined by the three moments 

52) 

of inertia (generally the moment of inertia tensor ). Moreover, when quantized spectra 
are considered, the tilted axis cranking gives a description for an isolated band just like 
the usual cranking does, but the wobbling motion as a whole corresponds to a multiple 
band structure, although we concentrated mainly on the first excited band in this paper. 

Some examples of the realistic calculations are also presented in this paper. The 
sign relation mentioned above leads characteristic zigzag behaviours in the E2 and Ml 
transition probabilities as functions of the rotational frequency, depending on equilibrium 
values of 7 deformation in the (e 2 , 7)-plane. This is quite analogous to the case of the 
B(M1) between the signature partner bands in odd- A nuclei. An important merit of the 
microscopic RPA formalism is that the three moments of inertia can be easily calculated, 
which is highly nontrivial. Their behaviour in three sectors of the (e2,7)-plane has been 
calculated in a rather simplified manner with neglecting selfconsistency of the mean-field 
potential. The results clearly show that their 7 dependence is neither irrotational nor 
rigid-body like. 

As a candidate of the nuclear wobbling motion which might be identified experimen- 
tally, we have also investigated the high-spin continuations of the so-called 7-vibrational 
bands with odd-spin in 182 0s, 124 Xe, and 126 Ba. The selfconsistent triaxiality deformation 
is somewhat different in each nucleus, and therefore the characteristic features of the col- 
lective vibrational modes are shown also to be different. Properties of the second-lowest 
RPA eigenmode, which is in some cases also rather collective but does not satisfy the 
criterion of the wobbling motion, are also discussed. 

Throughout this paper the Ml transition between the wobbling and yrast bands has 
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also been kept in scope. It should, however, be noticed that the expression of the B(M1) 
in the rotor model are not justified from the microscopic viewpoint. This is because the 
individual property of quasiparticle orbits are essential to understand the Ml transitions, 
which is not taken into account in the macroscopic model, in contrast to the case of the E2 
transitions, where the geometrical shape of the nucleus as a whole is mainly responsible. 
Although the Ml transitions from the quadrupole vibrational band are generally expected 
to be small at low-spin, it has been shown that this is not always the case at higher-spin 
region, where the rotation-aligned quasiparticles present in the vacuum. Note that the 
Ml transition is only possible if there exist the rotationally induced K mixing in the RPA 
phonon amplitudes. As a conspicuous example the high-spin continuation of the odd-spin 
7-vibrational band in 126 Ba is predicted to have much stronger Ml transition probability 
than the AI = ±1 E2, which is consistent with the observed branching ratio. This strong 
Ml transition, on the other hand, gives a possible reason why the AI = —2 transition that 
links above and below the band crossing in the odd-spin 7 band has not been observed. 

It is worthwhile to stress the importance of the collective vibrational motions of either 
the wobbling like or of the 7-vibrational character for the understanding the electromag- 
netic property of the odd- A nuclei. The characteristic property of the vibrational modes 
reflect strongly on the transition rates through the particle-vibration coupling effects, and 
therefore the study of the odd- A nuclei gives a valuable testing ground to clarify the char- 
acter change of the vibrational mode predicted in the present work. We hope that new 
generation of the large array of the crystal ball will provide us more detailed information 
of the electromagnetic transitions, which is necessary to confirm the predictions and to 
identify the nuclear wobbling motion if they exist. 
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Table 1. 

The first order coefficient of the rotational frequency dependence calculated 
by the RPA theory, a 7 in eq.(2.5) and the extracted intensity relation param- 
eters, Q t and q in eq.(2.9b). Here (g) and (7) attached to the theoretical q 
values mean that they are obtained by eq.(2.9b) using experimental moment of 
inertia of the ground band (J = J g ) and of the 7 band (J = J 7 ). 



Nucl. 


a-y/h [MeV- 1 ; 


Qt h) [eb] 


Q\ exp) [eb] 


q (th) {g) 


g(«0( 7 ) 


q(exp) 


162 Dy 


2.0 


0.37 


0.35 1 ) 


0.027 


0.024 




164 Dy 


1.9 


0.38 


0.34 2 ) 


0.023 


0.020 


0.021 2 ) 


164 Er 


2.0 


0.41 


0.38 3 ) 


0.031 


0.027 




166 Er 


1.8 


0.41 


0.42 4 ) 


0.024 


0.021 


0.022 4 ) 


168 Er 


1.3 


0.42 


0.37 5 ) 


0.017 


0.015 


O.OI8 5 ) 



1) Extracted from the B(E2 : 0+ -> 2+). ; 

2) Extracted from the x 2 fitting of the matrix elements M. (E2; AI = 0, —2) 1 ^ g 
(7 7 = 2, 4, 10) obtained by the Coulomb excitation experiment of Ref.. j4) 
There is an appreciable scattering of data from the relation, eq.(2.9a). The 
weights inversely propotional to the square of error bars of data are used 
for the x 2 fitting. 

3) Extracted from the B(E2 : 0+ -> 2+)f 5) 

4) FromRef.. 19) 

5) Extracted from the x 2 fitting of the matrix elements A4(E2; AI = 
0, ±1, ±2) 1 ^ g (Ij = 2, 3, 8) obtained by the Coulomb excitation experi- 
ment in Ref. ° 6) just in the same way as 2). 
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Figure Captions 



Fig.l The electric E2 transition amplitudes, eq.(2.4), in 164 Er microscopically calculated by 
the RPA formalism as functions of the rotational frequency. 

Fig. 2 Schematic figure depicting the analogy of the AI = ±1 E2 transitions between the 
yrast and the wobbling bands in even-even nuclei with the Ml transitions between 
the signature partner bands in odd nuclei. 

Fig. 3 Schematic figure depicting the relation between the triaxiality of the mean-field and 
the AI = ±1 E2 transitions from the wobbling band to the vacuum band. The (e2, 7) 
plane is classified into three regions, for which characteristic zigzag behaviours of 
B(E2) l \i=±i are shown in the middle panels. The transitions with stronger B(E2) , s 
are marked in the spectra (right panels). 

Fig. 4 The E2 transition amplitudes, Q y and Q z (left panel) and AI = ±1 B(E2) (right 
panel) for the lowest RPA solutions with signature a = 1 as functions of the rotational 
frequency in 182 0s. The calculational procedures are the same as those in §2 but the 
mean-field parameters are fixed as €2 = 0.21, and 7 = —18°, and A u = 0.74, A n = 0.98 
MeV, which are appropriate for the s-bands in 182 0s. 2j) As a reference the Weisskopf 
unit of B(E2) in this case is 0.0061 e 2 b 2 . 

Fig. 5 The effective moment of inertia as functions of the triaxial parameter 7 at hu) ro t = 0.25 
MeV. The results for -120° < 7 < -60° in 176 Hf, for -60° < 7 < 0° in 182 0s and for 
0° < 7 < 60° in 148 Gd are gathered in one panel. The mean- field parameters except 
7 are fixed as €2 = 0.27 and 0.22 for 176 Hf and 148 Gd, respectively, and A v = 0.6, 
A,,- = 0.6 MeV for both nuclei. The same parameters are used for 182 Os as Fig. 4. Note 
that 176 Hf and 148 Gd are known to be the nuclei with prolate and oblate non- collective 
rotations so that the selfconsistent 7 deformations in these nuclei are 7 = —120° and 
60°, respectively. The diabatic configurations of the vacuum states are chosen to be 
two-quasineutron four-quasiproton states for 176 Hf, two-quasineutron states (s-band) 
for 182 0s, and two-quasineutron and two-quasiproton states for 148 Gd. For 176 Hf and 
148 Gd, these configurations correspond to the high-X states with spin 20 + and 18 + ft, 
respectively, which are the yrast with (tt, a) = (+, 0) for both neutron and proton at 
their selfconsistent deformations. Unfortunately these high-X states do not coincide 
with any of the observed yrast isomers. 

Fig. 6 Routhians for the lowest lying state in the a = 1 sector as functions of the rotational 
frequency in 124 Xe. The calculated one is shown by the solid line and the experimental 
one, which is the routhian of the odd-spin members of the so-called 7 band, by the 
dashed line, and the observed yrast band by the short-dashed line. Experimental 
routhians are shown with respect to a reference band with J§ = 15.07i 2 /MeV, J\ = 
25.0^ 4 /MeV 3 , e s = 3.03 MeV and i s = 7.14%. Parameters used for the calculation 
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are e 2 =0.19, 7 = -45° and A„ = 0.9, A w = 1.1 MeV. 
Fig. 7 Calculated effective moments of inertia of the lowest lying state in the a = 1 sector 
in 124 Xe as functions of the rotational frequency. Parameters used are the same as in 
Fig.6. 

Fig.8 The E2 transition amplitudes, Q y and Q z (left panel) and AI = ±1 B(E2) (right 
panel) for the lowest RPA solutions with signature a = 1 as functions of the rotational 
frequency in 124 Xe. Parameters used are the same as in Fig.6. As a reference the 
Weisskopf unit of B(E2) in this case is 0.0037 e 2 b 2 . 

Fig. 9 The Ml transition amplitudes, iji y and /j, z (left panel) and B(M1) (right panel) for 
the lowest RPA solutions with signature a = 1 as functions of the rotational frequency 
in 124 Xe. Parameters used are the same as in Fig.6. As a reference the Weisskopf 
unit of B(M1) is 1.79 n 2 N . Note that the relative phase between the Ml amplitudes 
and the E2 amplitudes (Fig.8) is meaningful and gives the sign of the E2/M1 mixing 
ratio. 

Fig. 10 Excitation energies (left panel) and the E2 transition amplitudes (right panel) of the 
lowest and second-lowest states in the a = 1 sector in 124 Xe calculated at %uj rot = 0.3 
MeV as functions of 7. Parameters used for the calculation except 7 are the same as 
Fig.6. 

Fig. 11 The same as Fig.6 but for 126 Ba. Parameters used for the calculation are 62 =0.24, 
7 =0, and selfconsistently calculated pairing gaps are 1.16 > A„ > 1.03 MeV and 
0.76 > A w > 0.71 MeV for 0.4 < hu vot < 0.5 MeV. A reference band with J = 
14.0ft 2 /MeV, J x = 31.0ft 4 /MeV 3 , e s = 3.00 MeV and i s = 7.64ft is used here. Note 
that the 7-ray which links above and below the band crossing has not been observed 
(see the text). 

Fig. 12 The same as Fig.8 but for 126 Ba. Parameters used are the same as in Fig.ll. As a 

reference the Weisskopf unit of B(E2) in this case is 0.0038 e 2 b 2 . 
Fig. 13 The same as Fig. 9 but for 126 Ba. Parameters used are the same as in Fig.ll. 
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